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I.  INTRODUCTION 

The  intensive  work  performed  on  the  detailed  behavior  of  diodes  has 
brought  about  a  major  change  in  the  national  program  for  achievement  of 
Inertial  Confinement  Fusion  using  particle  beam  drivers.  This  program  has 
now  shifted  from  electron  beams  to  the  light  ion  beams  generated  in  diodes 
with  high  efficiency.  These  ions  will  be  focused  onto  medium  size  targets 
(about  1  cm  in  size)  from  many  ion  sources  in  order  to  achieve  the  needed 
high  power  (about  100  terrawatt)  for  10  to  20  ns.  The  detailed  understanding 
of  how  such  diodes  turn  on  ion  sources  in  order  to  generate  about  1017  ions 
from  a  single  diode  is  of  crucial  Importance.  The  time  dependence  of  the 
ion  current,  the  type  of  ion  species,  the  magnetic  field  distribution,  anode 
and  cathode  plasma  motion  (affecting  ion  trajectories  and  thus  focussing) 
have  all  been  investigated.  The  effects  of  different  diode  materials, 
geometries,  etc.  are  under  continuing  study.  A  combined  theoretical  and 
experimental  research  program  of  the  physics  of  light  ion  beams  is  presently 
supported  by  JAYCOR.  Many  of  the  results  obtained  at  the  Naval  Research 
Laboratory,  Sandla  National  Labs,  Physics  International  and  Cornell 
University  have  been  explained  successfully  by  the  diode  theory  developed 
by  JAYCOR  personnel.  The  physics  of  electron  flow  in  very  large  diodes 
(of  the  order  of  1  meter)  as  well  as  of  the  anode  behavior  (such  as  anode 
plasma  formation)  Is  of  crucial  importance  to  the  question  of  power  transfer 
to  small  area  loads  (such  as  small  diodes  or  plasma  loads). 

In  diode  geometries  that  converge  to  small  radii,  there  can  be  signi¬ 
ficant  electron  loss  to  the  anode  in  the  vacuum  feed  which  may  turn  on  an 
anode  plasma  there.  In  addition,  the  presence  of  this  ion  source  is  another 
cause  for  energy  drain  on  the  power  flow  since  the  ion  current  is  enhanced 
by  a  large  factor  due  to  electron  beam  near  the  anode. 


The  three  subjects,  small  size  diodes  (1  mm)  in  order  to  generate  very 
large  ion  current  density,  medium  size  diodes  (10  cm)  for  efficient  ion  beam 
generation  via  the  Pinch  Reflex  Diode. and  large  diodes  (-1  m)  were  the 
principle  areas  Investigated  and  are  concentrated  on  in  this  report.  In 
particular  this  report  explains  the  physics  of  the  backscatter  diode  and  how 
a  runaway  phenomenon  can  generate  very  large  current  density  Ion  beams  which 
could  be  used  in  beam  target  studies.  There  follows  an  examination  of  the 
physics  of  how  anode  plasmas  are  generated  using  various  materials  and  the 
question  of  electron  flow  in  a  pinch- reflex  diode  (PRD).  The  report  ends  with 
a  review  of  the  role  of  electron  and  ion  flows  in  large  area  diodes  driving 
plasma  loads v 


II.  SUMMARY  OF  WORK 

The  work  performed  on  the  physics  of  electron  and  ion  flows  in  diodes 
was  a  combination  of  an  analytical  study  and  complimentary  numerical  simulation 
using  Particle  In  Cell  codes.  These  codes,  developed  during  the  last  few 
years,  have  been  checked  by  comparing  their  predictions  to  the  experimentally 
observed  variables.  They  were  found  to  be  in  close  agreement.  The  theo¬ 
retical  analysis  was  based  on  general  descriptions  of  relativistic  electron 
drift  motions,  anode  material  properties,  plasma  generation  in  different 
gases  and  Monte-Carlo  results  for  electron  scattering  in  materials. 

In  the  area  of  the  backscatter  diode  it  was  discovered  that  even  a 
single  backscatter  (per  electron)  with  very  high  energy  loss  (more  than  half 
of  the  diode  voltage  lost  during  the  backscatter  process)  results  in  no 
steady  state  operation.  This  explains  phenomena  observed  in  a  few  labora¬ 
tories  studying  diodes  immersed  in  very  large  axial  magnetic  fields.  The 
backscatter  diode  is  contrasted  with  the  reflex  diode  and  the  similarities 
of  their  load  characteristics  for  the  external  generator  are  pointed  out. 

We  also  show  that  very  high  current  densities  (>107A/cm2)  can  be  achieved 
using  large  (>105G)  axial  magnetic  fields. 

In  the  area  of  the  Pinch-Reflex-Diode  the  time  scale  for  anode  plasma 
generation  using  different  anode  materials  was  investigated.  The  case  of 
pinch  formation  without  ion  flow  was  investigated  using  the  PIC  code.  The 
fraction  of  electrons  in  the  pinch,  diode  impedance  and  energy  deposition 
in  the  foil  were  computed.  The  ion  current  was  then  turned  on  In  the  pinch 
area  and  out  to  larger  radii  and  was  found  to  scale  linearly  with  r  (that  is 
J-f(r)  —  r“^).  Based  on  the  deposition  profile,  turn-on  times  of  anode  plasma 
were  computed  and  compared  to  experimental  results.  The  Impedance  of  the 


diode  was  then  computed  with  ion  current  flowing,  showing  a  decrease  by  a 
factor  of  1.5  relative  to  the  same  diode  operating  before  ion  flow  exists. 

On  the  topic  of  large  diodes  and  imploding  plasmas,  the  role  of  "mini¬ 
diodes"  created  in  the  plasma  was  explained.  These  psuedo-diodes  are 
responsible  for  carrying  a  large  fraction  of  the  current  across  the  vacuum 
gaps  in  sausage  instability  regions.  The  dynamics  of  implosion  is  followed 
from  early  plasma  formation  on  wire  arrays  to  the  driving  forces  that  col¬ 
lapse  the  plasma.  The  stability  of  the  run  in  phase  is  attributed  to  the 
low  inductance,  return  current  paths. 

The  relativistic  electrons  that  flow  into  the  plasma  load  from  the 
vacuum  feed  lose  a  significant  fraction  of  their  energy  to  the  plasma 
electrons  due  to  orbit  lengthening  caused  by  the  magnetic  field  there. 

When  currents  of  a  few  megamperes  flow  in  the  plasma,  the  resulting  self 
field  enhances  the  stopping  power  by  more  than  two  orders  of  magnitude  over 
the  classical  value. 

The  very  large  diode  (-1  m)  that  feeds  power  to  the  plasma  load  was 
studied  regarding  the  problem  of  magnetic  insulation.  In  particular,  an 
estimate  was  obtained  for  the  fraction  of  diode  power  lost  to  electron  flow 
in  vacuum.  The  effect  that  the  electron  launch  geometry  at  very  large  radii 
as  well  as  the  Impedance  of  the  generator  have  on  that  flow  was  pointed  out. 
The  high  impedance,  high  power  diodes  prove  to  be  the  most  susceptable  to 
power  losses  from  either  electron  or  ion  flows  in  vacuum.  Ion  flows  are 
computed  for  typical  -in  Impedance  diodes.  The  time  scale  for  generating 
these  flows  is  then  pointed  out.  If  the  present  simple  designs  were  scaled 
to  higher  powers  {to  say,  30  TW)  more  than  half  the  power  could  be  lost. 
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III.  THEORY  OF  FLOWS  IN  DIODES 

In  this  section,  three  separate  diodes  are  treated:  the  backscatter 
diode,  the  PRO  and  the  diode  feed  for  plasma  loads.  These  diodes  can  be 
conveniently  characterized  by  their  sizes  which  are  roughly  1  mm,  10  cm  and 
1  meter  respectively.  These  sizes  are  dictated  by  the  applications  of  these 
diodes.  The  backscatter  diode  may  be  a  tool  to  generate  currents  of  only 
105A  but  from  small  (10  cm2)  areas  so  that  very  high  current  density 
fluxes  may  be  obtained  and  used  in  beam  target  studies  at  current  densities 
that  exceed  by  an  order  of  magnitude  the  present  10®  A/cm2.  The  PRD  is  a 
diode  for  generating  ion  currents  very  efficiently  (60-70%)  at  high  currents 
of  10®A  and  directing  these  ions  to  a  focus  with  a  small  launching  angle 
(half  cone  of  10°).  The  size  of  such  a  diode  cannot  be  of  the  order  of  1 

cm  since  then  its  impedance  will  be  much  larger  than  the  desired  in  (which 

would  require  a  gap  of  0.5  mm  which  closes  too  quickly  by  plasma  motion 
relative  to  the  needed  ~100  ns  pulse).  It  also  cannot  be  excessively  large 
(say  30  cm  radius)  sirce  then  focussing  distances  of  150  cm  are  required  in 

assuming  a  10°  launching  angle.  At  such  distances  a  focal  spot  of  1  cm 

requires  ion  beam  quality  of  1/3°.  Since  present  beams  have  divergences  of 
1.5°  -  3°  the  diode  diameter  has  to  be  about  10  cm.  The  size  of  the  diode 
feeding  the  plasma  load  is  constrained  by  impedance  considerations  for 
the  generator  and  by  water  breakdown  requirements  for  a  given  power  of  a 
few  terrawatts.  These  constraints  lead  to  one  meter  sizes. 


Discussions  of  the  details  for  each  case  follow. 


1.  Backs cat ter  Diode. 


When  electrons  emitted  from  a  cathode  strike  a  high  atomic  number 
anode, a  large  fraction  of  the  electrons  are  backscattered  into  the  diode. 
During  the  scattering  in  the  material  they  loose  energy  and  emerge  back  into 
the  diode  with  an  energy  eV  <  eV  where  V  is  the  diode  voltage  and  e  the 

5  V  w 

electron  charge.  In  addition,  the  momentum  distribution  of  the  scattered 

electron  fills  nearly  a  hemisphere.  In  the  presence  of  very  large  axial 

magnetic  fields  the  perpendicular  momentum  is  conserved  during  electron 

motion  in  the  diode  and  the  parallel  momentum  serves  to  drive  the  electrons 

back  towards  the  cathode.  The  point  of  closest  approach  is  reached  when 

0.  Given  and  P„  at  the  anode  for  a  backscattered  electron  it  will 

★ 

go  back  to  a  point  having  a  voltage  V  (the  cathode  is  at  V  »  0,  the  anode 
at  V  *  VQ).  The  electrons  leave  the  anode  with  kinetic  energy  eVs.  This 
implies 

eV. 

y( anode)  =  1  +  |  g-|2  | 

which  is  related  to  their  momentum  by  the  relativistic  relation 


Y 


2 


1  + 


+  p,; 


As  the  electrons  flow  back  into  the  diode  from  VQ  at  the  anode,  they 

loose  kinetic  energy  and  gain  potential  energy;  for  the  relativistic  case 

they  loose  only  P„  (in  fact,  vA  is  not  conserved  but  rather  increases!), 

* 

and  at  any  voltage  V  <  V  <  VQ  the  relativistic  factor  is 


from  which  one  may  find  P„(at  V).  The  value  of  V  is  given  by  setting  ?H 
★ 

(at  V  )=  0  in  the  momentum  equation 


(’  +  *VT!) 


1  + 


e(V$  -  VQ  +  V  ) 


mQc2 


For  example,  if  the  electron  lost  no  energy  in  the  scattering  process 
V$  =  V0  and  if  it  got  scattered  with  P£  =  P*  =  P*  then  P2  -  2  P,*  =  ^  P* 
where  P0  is  the  total  momentum.  The  value  of  PQ  is  related  to  VQ  by 


eV, 


Po 


(1  +  Ijjf^rl)2  =  1  +  •  After  some  algebra  one  finds  V 


eV 

me2 

o 


=  1  - 


1  + 


giving  about  0.35  MV  for  VQ  *  0.5  MV.  The  typical  electron  does  not  approach 
the  cathode  but  stays  near  the  anode  side  (that  is  its  voltage  is  hot  less 
than  0.25  MV).  The  first  problem  to  realize  is  that  the  electron  accumula¬ 
tion  in  the  diode  before  ions  are  generated  causes  a  reduction  in  the 
electric  field  near  the  cathode  and  suppression  of  emission.  The  current 
produced  by  N  backscatters  can  be  computed  approximately  by  using  each 
stream  of  reflected  electrons  as  having  an  amount  of  energy  e(VQ  -  AV$) 
after  the  S-th  reflection  (see  Fig.  1-1).  A  solution  is  sought  for  the 
current  density  emitted,  dg.  In  a  diode  gap,  D,  with  diode  voltage  VQ. 

For  simplicity,  relativistic  effects  are  disregarded.  Poisson's  equation 
in  the  axial  direction  -  x  is: 


where  s  a  0  Is  the  first  stream  of  electrons  from  the  cathode*  s  a  1  is  the 
first  backscattered  flow,  etc.  The  density  Is  double  for  s  >_  1  since  that 

e+1  c 

stream  returns  symmetrically  to  the  anode.  Notice  that  0g  <.  since  we 
assume  that  electrons  can  only  be  lost  (realistically,  secondary  electrons 
are  also  generated  and  for  special  cases  an  electron  flow  back  after  scatter¬ 
ing  could  be  enhanced). 

The  velocity  Is  given  by 


\  m(v*)2  a  e(<J>  -  AV$) 


The  Poisson  equation  can  now  be  written  in  simplified  form 


d2<t>  _  „o,o  v  Ks  , 
d#  K  Je 

which  is  of  the  form  =  F(<J>)  that  can  be  integrated  once  by  multiplying 
by  Hx  t0  ($<)  *  gives  by  integration 


?<3f>2 -r<f>!l  ■/«♦>«♦  . 

0 

For  the  case  of  space  charge  limited  emission  from  the  cathode  with  a  cathode 

temperature  much  less  than  the  diode  voltage  (typically  Tcathode  <_  10  eV,  while 

V  «~1  MeV)  and  with  a  small  work  function  for  the  plasma,  the  electric  field 
o 

In  the  plasma  is  a  few  kv/cm  compared  to  the  ~1  Mv/cm  In  the  diode. 

Matching  the  field  at  the  plasma/electron  beam  interface  on  the  cathode  we 

may  thus  set  ^  |  *  0.  The  equation  for  ^  may  now  be  Integrated  again 
xao 
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and  setting  the  limits  on  the  anode, $  *  VQ  and  x  *  D,  one  finds  the  single  un 

known  Notice  that  the  integration  on  the  sum  of  flows  starts  for  the 

s-th  term  at  $  *  AV$  (since  no  electrons  flow  closer  to  the  cathode  than 

AV$)  thus  no  Imaginary  values  are  included  in  our  treatment.  The  result 

after  straightforward  but  cumbersome  integration  is  that  the  current  density 

is  given  by  the  Child-Langmuir  expression  with  a  suppression  factor  that 

AV  s+1 

depends  on  the  values  of  tt-=-  and  =- —  .  This  treatment  will  not  be  given 

vo  Js 

here  since  It  adds  no  special  understanding  of  the  physics  of  the  diode. 

s+1  s  AVs  s 

An  approximate  form  was  derived  for  the  case  where  J  /J  «  1  and  y-2-  «  £ 

o 

which  is  for  N  »  1: 

2.3  V03/2(MW  ^  w 
e  D2  Not5* 

The  current  density  is  thus  reduced  from  the  no  scattering  case  and  for 
N  *  10  It  is  about  4  times  lower.  The  diode  perveance  is  thus  lower  if 
one  uses  high  Z  material  anodes  Instead  of  low  Z.  Initiation  of  the  ion  flow 
from  the  anode  side  adds  another  term  in  the  Poisson  equation  for  the  ion 
space  charge  and  gives  us  one  more  boundary  condition  for  the  emission  law. 

If  anode  plasma  Is  assumed  then  for  similar  reasons  to  those  on  the 
cathode  side  the  electric  field  is  zero  on  the  anode.  There  are  thus  two 
unknown  current  densities,  Jfi  and  emlnatlng  from  the  cathode  and  anode 
respectively.  In  addition,  electrons  are  being  backscattered  from  the 
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anode.  Instead  of  solving  the  general  case  with  its  many  unessential 
details,  a  single  backscatter  of  all  electrons  is  assumed  with  a  loss  of 
energy  AV.j.  Notice  that  the  ions  (assumed  singly  charged  with  mass  M)  move 
from  the  anode  and  thus  their  velocity  is  given  by 

^MV*  -  e  (VQ  -  V)  . 

In  the  region  between  $  *  AV-|  and  $  s  Vfl,  the  space  charge  is  composed  of 

”  ^e 

three  flows:  the  primary  electron  flow  with  density  (2  e'/mfe  *  the  back* 

-2  J 

scattered  flow  with  density  '^'e/n^V'^-AV  )**  *  an<*  ^on  ^ow  a 

1 

[i.  e/M)^"(V  -4)]%  *  Poisson's  equation  may  thus  be  written  as 


4-rrJ 


(2e/m) 


7*  +  (<f>  -  AV-j  )**  ’  (VQ  -  ♦)**' 


(VV  (M/m)1* 


Integrating  the  equation  once  from  x=o  by  the  same  technique  used  earlier 

and  using  *  0  one  finds: 

ax  x*o 

I  •  TT^  )2  *  .«♦  -  ‘V*  +  2[‘vo  -  ♦>'  - 

*cji/je)  (w1* 

The  L.H.S.  is  also  zero  at  x*D  where  4>*V0  giving  an  equation  for  Jj/Je: 

V*  +  2 (VQ  -  AV-|  )**  -  (Ji/Je)  (M/m)1* 
or: 

k  k 

(M/m)**  J./J  -1+2(1-  j+r 

0 

which  shows  enhancement  over  the  bipolar  ratio  by  a  factor  equal  to  the  R.H.S. 

11 


The  above  expression  may  now  be  inserted  into  the  equation  for 

and  one  may  proceed  to  solve  by  another  integration  for  <p(x) .  Then,  since 

♦(D)  "  V-,  on®  can  find  the  last  unknown:  The  surprise  is,  however, 

0  AV,  e  ^ 

that  depending  on  there  may  or  may  not  be  a  real  solution  to  ^  through 

Its  equation: 

?  <$>'  ■  ^7^  |2  *H  *  «<♦  -  * 2  K  -  ♦>'s  -  O 

*[1  +  20  -?>"]}  • 


It  is  interesting  to  note  how  the  R.H.S.  behaves  as  <p  goes  from  0  to  VQ. 
It  starts  from  0  and  Increases  until  at  ♦  *  AVj  it  equals 


(AV,)**  +  [(VQ  -  AV-j )**  -  V*]  [l  +  2(1  -  JV] 


2F 


AV, 

The  quantity,  2F,  can  be  expanded  for  n-2-  *  e  «  1  and  for 

vo 

to  yield 


5  «  1 


IV^e*  (1  -  j  e*5)  for  e  «  1 

61*  (j  6**”-  1 ) .  for  5  «  1  . 

This  gives  F  >  0  for  c  but  F  <  0  for  6. 

It  can  thus  be  concluded  that  very  small  energy  losses  (s  case) 


give  steady  state  solutions  but  very  large  losses  (6  case)  end  up  with  no 
solution  of  steady  electron  and  Ion  flow.  One  would  then  like  to  know  where 
the  dividing  line  is  and  what  would  happen  physically  when  the  steady  state 


goes  out  of  existence?  The  first  question  is  easy  to  answer:  set  F  *  0 
and  solve  for  AV-j  to  find  that  AV^  »  j  VQ.  Since  it  has  already  been 
shown  that  typical  very  high  z  backscatter  gives  AV^  >  ^  VQ,  it  can  be 
concluded  that  no  steady-state  solutions  exist  for  these  materials.  The 
time-dependent  solution  for  general  materials  is  beyond  the  scope  of  this 
work  but  it  should  be  realized  that  low  z-material  backing  a  thin  layer  of 
high  z  can  act  as  an  absorber  of  electrons  and  thus  bring  back  a  steady 
state  solution.  The  second  question  can  be  treated  only  heurlstically. 

The  ion  space  charge  Is  enhanced  by  electrons  present  near  the  anode.  It 
flows  towards  the  cathode  where  it  finds  only  the  primary  electron  flow 
without  backscattered  electrons  (.$  <  AV-|).  This  ion  space  charge  is  too 
large  for  the  existing  electron  charge  present  and  it  thus  causes  increased 
emission  from  the  cathode.  After  the  newly  emitted  electrons  arrive  and 
backscatter  off  the  anode  they  generate  a  larger  ion  space  charge  which 
in  turn  after  a  transmit  time  to  the  cathode  repeats  the  process  causing 
the  total  current  to  runaway.  The  current  density  is  limited  only  by  self 
fields  that  change  the  electron  flow.  Thus  very  large  axial  magnetic  fields 
are  needed  for  small  area  cathodes.  The  total  current  is  limited  by  the 
generator  and  the  diode  voltage  drops  by  inductive  effects.  When  the 
voltage  is  very  low  only  a  few  electrons  backscatter  and  the  diode  can  find 
a  steady  state  (saturation  of  runaway).  The  current  may  stop  increasing 
and  this  inductive  kick  will  bring  the  voltage  back  to  the  threshold  between 
existence  and  non-existence  of  steady  state  in  the  electron/ion  flow  problem. 
This  process  is  very  similar  to  the  reflex  diode  but  no  double  diode 
arrangement  is  needed.  The  mathematics  are  very  similar  and  the  processes 
observed  are  almost  the  same.  They  are  summarized  in  Figure  1-2  where  the 
time  dependence  of  current  and  voltage  are  shown  as  a  function  of  time. 
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Figure  1-2 


RANGE  IN  ANODE  FOIL  (NEED  ~FOUR 
APPEARANCES) 

NO  AREA  OR  A-K  GAP  CONSTRAINT 
ON  TOTAL  DIODE  CURRENT 


2.  Pinch-Reflex  Diode. 

a.  Anode  Plasma  Generation.  The  problem  of  anode  plasma  genera¬ 
tion  using  different  anode  materials  Is  fairly  complicated  when  absorbed 
gases  are  Included  In  the  process  of  plasma  generation.  For  clean  surfaces 
of  materials  such  as  copper,  aluminum  and  polyethylene,  first  the  simple 
sol Id-to- liquid  then  the  liquid-to-gas  phase  transition  and  finally  the 
gas  Ionization  Into  plasma  have  to  be  treated.  In  the  first  heating  phase, 
relativistic  electrons  deposit  their  energy  into  the  anode  material  via  the 
classical  single  particle  slowing  down  process.  The  typical  stopping  power 
of  about  1.5  MeV/g/cm2  Is  enhanced  by  electron  scattering  in  thick  materials. 
This  scattering  effect  was  a1rea4y  treated  by  Monte-Carlo  methods  and  the 
enhancement  factors  are  given  In  the  NBS  manuscript  by  L.V.  Spencer  on 
Energy  Dissipation  by  Fast  Electrons  (the  results  for  the  1  MeV  typical  for 
Ion  diodes  is  used  herein).  For  polyethylene  the  stopping  power  is  similar 
to  that  of  carbon  which  is  1.6  MeV/g/cm2.  This  is  enhanced  by  scattering 
by  a  factor  of  1.28  resulting  in  2  MeV/g/cm2  on  the  surface  of  the  material. 
For  aluminum,  1.51  MeV/g/cm2  is  enhanced  by  a  factor  of  1.87  giving  about 
2.8  MeV/g/cm2.  For  copper  1.29  MeV/g/cm2  is  enhanced  by  3.64  giving  4.7 
MeV/g/cm2.  All  of  the  above  are  for  a  flow  of  electrons  perpendicular  to 
the  anode  surface  with  all  the  backscattered  electrons  returning  to  the 
anode  (an  assumption  justified  by  the  large  electric  field  in  the  diode 
where  anode  heating  takes  place).  The  next  question  is  how  much  charge 
per  unit  area  Is  needed  In  order  to  bring  anyone  of  the  above  materials 
from  room  temperature  to  evaporation?  Assume  that  an  energy,  E  (specified 
In  joules/gram).  Is  needed  for  evaporation.  The  energy  deposited  during 
t(seconds)  by  a  current  density  J(A/an2)  impinging  on  the  surface  of  a 
material  with  stopping  power  SCeV/g/cm2)  Is  given  by  the  product  JtS  (joules/ 
gram).  In  order  to  reach  the  value  of  E  for  a  given  material  one  has  to  have 
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J  t  -  E/S  . 


Therefore,  It  is  the  charge  delivered  per  unit  area  by  the  1  MeV  electron 
beam  that  determines  whether  or  not  the  material  is  evaporated.  The  values 
for  E  are  given  by  Hultgren  R. ,  Orr  R.,  Anderson,  P.D.,  and  Kelly,  K.  K. 
in  Selected  Values  of  Thermodynamic  Properties  of  Metals  and  Alloys  (New 
York,  John  Wiley  and  Sons,  Inc.,  1968).  The  values  for  aluminum  and  copper 
are  given  as  14  and  6  kJ/g  respectively.  A  separate  calculation  is  needed 
for  polyethylene.  After  evaporation  the  molecular  bonds  must  be  broken 
before  plasmas  of  hydrogen  and  carbon  can  be  generated  from  the  (CH2)n 
polymers.  The  reaction  CH2  +  e  -*■  CH2  +  2e  requires  11.9  eV  or  80  kJ/g. 
Similar  values  are  found  for  CH2  +  e  -*•  CH  +  H+  +  2e.  The  stripping  of 
electrons  from  these  molecules  is  very  similar  to  that  from  water  (12.6 
eV  for  H20  +  e  H20+  +  2e).  This  information  can  be  found  in  Herzberg, 
Infrared  and  Raman  Spectra  (page  528).  Using  the  above  values  we  find  that 
J  t  Is  equal  to  4x10  ,  5x10  ,  1.3x10  Coulomb/cm2  for  polyethylene, 

aluminimum  and  copper  respectively.  For  a  time  scale  of  say  10  ns  the 
needed  current  densities  for  evaporation  are  4x10s,  5x10s,  1.3x10s  A/cm2 
for  the  above  materials.  These  results  explain  why  It  is  that  conducting 
anodes  coated  with  (CH2)h  or  CD2)n  do  not  generate  Ion  currents  until  very 
late  into  the  pulse  (about  50  ns  assuming  0.8  x  10s  A/cm2).  The  aluminum 
takes  an  order  of  magnitude  less  time  If  it  was  first  baked  to  remove  the 
absorbed  gases  (a  technique  used  in  experiments  at  Sandla  Labs,  Albuquerque 
and  at  Physics  International),  Without  prebaking  the  gas  desorption  process 
occurs  at  a  low  temperature  early  In  the  pulse  and  the  expanding  gas  is 
Ionized  by  an  avalanche  process  triggered  by  the  diode  electric  field.  This 
Ionization  relies  on  the  100  eV  electrons  generated  during  the  avalanche 
rather  than  the  1  MeV  relativistic  electron  beam.  The  cross  section  for 
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ionization  by  a  100  eV  electron  is  three  orders  of  magnitude  larger  than 
that  by  a  1  MeV  electron  and  thus  currents  of  10“  A/cm2  can  generate  the 
plasma  in  nanosecond  time  scales. 

The  heating  process  for  a  pinch  reflex  diode  is  enhanced  by  the 
electron  reflexing  through  the  foil.  The  performance  of  a  given  diode  is 
initially  investigated  assuming  pure  electron  flow  in  the  diode.  The 
steady-state  electron  flow  pattern  from  the  cathode  to  the  anode  along 
with  the  self-consistent  fields  are  modeled  by  means  of  a  Particle  In  Cell 
(PIC)  code  DI0DE2D  that  has  been  developed  over  the  past  few  years.  This 
code  was  used  to  simulate  the  basic  Pinch  Reflex  Diode  (PRD)  geometry  shown 
in  Figure  2-1.  The  electrons  reflex  through  the  anode  foil  due  to  the 
azimuthal  magnetic  field  created  by  the  current  flow  through  the  button 
connecting  the  thin  foil.  They  are  again  ejected  from  the  anode-cathode  gap 
by  the  diode  electric  field.  A  typical  electron  flow  pattern  is  shown  in 
Figure  2-2.  The  electrons  flow  to  the  small  central  conductor  (0.5  cm 
radius)  starting  from  the  hollow  cathode  (6.4  cm  radius)  with  hollow  shank 
thickness  of  0.7  cm  projecting  axially  out  from  a  solid  cathode.  The 
cathode  separation  from  the  anode  foil  was  0.32  cm  (giving  a  diode  aspect 
ratio,  R/D  »  20).  The  gap  between  the  foil  and  the  solid  anode  was  also 
0.32  cm.  In  addition,  the  length  of  the  hollow  cathode  shank  projecting  from 
the  solid  cathode  was  changed  from  zero  to  0.32  cm  and  no  significant  effects 
on  electron  flow  were  observed,  Since  the  electric  equipotentlal  lines  are 
compressed  near  this  anode  surface  It  Is  not  surprising  that  the  exact 
hollow  cathode  shank  length  projecting  Into  the  A-K  gap  has  little  effect 
on  the  E  field.  In  addition,  two  separate  diode  voltages  were  studied,  750 
and  1500  kV.  For  750  kV  the  total  electron  current  was  335  kA  (2.2  51).  The 
total  current  with  ions  and  electrons  was  540  kA  (1.4  0),  a  significant  drop 
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PINCH-REFLEX  DIODE  SIMULATION 


in  impedance.  This  drop  occurs  because  the  currents  are  self  magnetically 
limited  (unlike  the  one  dimensional  case  where  a  factor  of  two  drop  in 
impedance  would  result  due  to  ion  space  charge).  For  the  1500  kV,  the 
results  show  the  same  Impedances.  The  fraction  of  electron  current  that 
pinched  on  the  central  anode  conductor  was  about  40%  at  these  impedances. 

When  the  impedance  was  increased  to  about  20  £1  (at  7  MV)  the  fraction  in  the 
pinch  was  nearly  100%.  There  exists  a  non-continuous  transition  between 
very  strong  and  efficient  vacuum  pinching  and  the  pinch  reflex  behavior. 

It  is  suspected  that  electrons  are  allowed  to  migrate  via  drifts  in  the 
low  Impedance  cases  while  they  execute  figure  eight  orbits  Cthe  nonadiabatic 
case)  at  the  high  impedance.  The  Alfven  current  forms  a  crude  division  be¬ 
tween  these  two  regimes.  The  Alfven  solution  does  not  include  electric 
fields  nor  does  it  solve  fields  self  consistently  in  the  diode.  Thus  no 
theory  presently  exists  to  characterize  this  type  of  diode  behavior. 

3.  Imploding  Plasmas  and  Diodes. 

A.  General . 

The  contributions  to  the  physics  of  imploding  plasmas  can  be 
summarized  in  three  areas: 

0)  Dynamics  of  implosion. 

C2)  Relativistic  electron  beams  in  the  plasma. 

(.3)  Problems  in  magnetic  insulation  of  radially  converging  feeds. 

In  the  area  of  the  dynamics  of  Implosion  we  point  out  (a)  the  importance 
of  early  time  expansion  of  the  wire  plasmas,  (b)  a  possible  mechanism  that 
acts  very  early  and  (c)  how  to  affect  it  by  the  polarity  of  the  wire  holder 
geometry.  The  actual  motion  of  the  wires  is  shown  to  be  dictated  at  early 
time  by  the  repulsive  forces  from  the  return  current  rods  and  at  late  time  by 


the  attracting  forces  of  the  other  wires.  The  stability  question  of  the 
run-in  phase  is  thus  divided  into  two  parts;  the  repulsive  rods  play  an 
important  role  in  early  time  and  a  specific  geometry  is  suggested  for  stable 
run  in. 


The  behavior  of  relativistic  electrons  in  the  imploding  plasma  is 
treated  in  order  to  find  their  energy  losses  inside  the  plasma  compared  to 
energy  losses  on  the  anode.  The  magnetic  field  forces  very  tight  long  orbits 
to  the  relativistic  electrons  forcing  more  than  two  orders  of  magnitude 
increase  of  the  stopping  power  of  the  electrons  in  the  plasma.  Typical  cases 
are  treated  showing  that  present  experimental  arrangements  allow  most  of  the 
energy  in  the  electron  beam  to  be  absorbed  in  the  plasma.  About  one  Alfven 
current  is  allowed  to  be  deposited  in  the  anode  in  agreement  with  most 
experimental  data. 

The  problems  associated  with  magnetic  insulation  problems  in  the  feed 
of  radially  converging  structures  show  very  different  behavior  for  negative 
vs.  positive  polarity  in  the  inner  electrode.  Major  losses  associated  with 
negative  polarity  are  reduced  in  positive  polarity  operation.  For  very  high 
power  (30  TW)  the  example  treated  shows  that  new  physics,  not  observed  at  the 
5TW  level,  will  dominate  the  particle  flows  in  the  feed  reducing  the  fraction 
of  current  in  the  plasma  by  about  a  factor  of  2. 

B.  Dynamics  of  Implosion. 

The  dynamics  of  wire  arrays  driven  by  multi terawatt  generators  was 
Investigated  from  a  new  point  of  view  regarding:  1)  prevention  of  early  wire 
expansion,  2}  radial  collapse  of  the  wire  array  and  3)  stability  of  the  run 
In  wire  array.  The  three  areas  are  connected  to  each  other  since  the  early 
wire  dynamics  dictates  the  plasma  density  formed  which,  in  turn,  is  highly 


the  attracting  forces  of  the  other  wires.  The  stability  question  of  the 
run-in  phase  is  thus  divided  into  two  parts;  the  repulsive  rods  play  an 
important  role  is  early  time  and  a  specific  geometry  is  suggested  for  stable 
run  in. 

The  behavior  of  relativistic  electrons  in  the  imploding  plasma  is 
treated  in  order  to  find  their  energy  losses  inside  the  plasma  compared  to 
energy  losses  on  the  anode.  The  magnetic  field  forces  very  tight  long  orbits 
to  the  relativistic  electrons  forcing  more  than  two  orders  of  magnitude 
increase  of  the  stopping  power  of  the  electrons  in  the  plasma.  Typical  cases 
are  treated  showing  that  present  experimental  arrangements  allow  most  of  the 
energy  in  the  electron  beam  to  be  absorbed  in  the  plasma.  About  one  Alfven 
current  is  allowed  to  be  deposited  in  the  anode  in  agreement  with  most 
experimental  data. 

The  problems  associated  with  magnetic  insulation  problems  in  the  feed 
of  radially  converging  structures  show  very  different  behavior  for  negative 
vs.  positive  polarity  in  the  inner  electrode.  Major  losses  associated  with 
negative  polarity  are  reduced  in  positive  polarity  operation.  For  very  high 
power  (.30  TW)  the  example  treated  shows  that  new  physics,  not  observed  at  the 
5TW  level,  will  dominate  the  particle  flows  in  the  feed  reducing  the  fraction 
of  current  in  the  plasma  by  about  a  factor  of  2. 

B.  Dynamics  of  Implosion. 

The  dynamics  of  wire  arrays  driven  by  multi terawatt  generators  was 
investigated  from  a  new  point  of  view  regarding:  1)  prevention  of  early  wire 
expansion,  2}  radial  collapse  of  the  wire  array  and  3)  stability  of  the  run 
In  wire  array.  The  three  areas  are  connected  to  each  other  since  the  early 
wire  dynamics  dictates  the  plasma  density  formed  which,  in  turn,  is  highly 
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Incompressible  when  assembled  on  axis.  This  density  is  the  prime  factor 
when  it  comes  to  consider  radiation  from  the  plasma.  In  this  work,  no 
reference  is  made  to  quantitative  radiation  but  it  is  assumed  that  high  plasma 
densities  are  desirable  for  specific  outputs.  This  work  will  thus  concentrate 
on  effects  that  cause  low  plasma  density,  what  is  their  origin  and  how  to 
avoid  them.  Included  are  specific  recommendations  for  experiments  to  be 
performed  on  the  large  generators. 

(1)  Early  Wire  Expansion. 

In  previous  work  it  was  shown  that  the  electric  field  enhance¬ 
ment  on  the  small  radius  wires  in  the  presence  of  MV  voltages  causes  field 
emission  currents  to  dominate  over  currents  running  inside  the  solid  wire. 

This  causes  fast  wire  heating  and  expansion.  Later  in  time  magnetic  fields 
outside  the  wire  plasma  grow  (during  current  rise)  and  start  to  compress  the 
plasma.  Skin  current  models  show  wire  plasma  compression.  The  model  contains 
a  very  dominant  assumption  of  current  in  the  skin  only.  In  reality,  magnetic 
field  diffuses  in  because  of  small  wire  plasma  radii  (<0.5  mm)  and  finite 
resistivity  early  In  time  (when  plasma  temperatures  are  ~10  eV).  This 
computation  Is  shown  in  Appendix  A.  The  main  effect  of  the  magnetic  field 
that  diffused  Inside  (say  at  half  the  peak  current)  is  that  when  the  plasma 
gets  hotter  It  freezes  inside  the  plasma  and  during  plasma  compression  it 
Increases  rapidly  acting  as  a  pressure  against  compression.  This  is,  in 
particular,  the  case  as  the  wire  plasmas  run  radially  Inward  to  assemble  on 
axis.  During  this  run-in  the  total  diode  Inductance  Increases  very  fast  and 
the  currents  through  the  multiple  wires  drop  (part  of  this  current  is  then 
carried  outside  the  wire  array)  thus  reducing  the  pressure  that  holds  each 
wire  plasma.  During  this  time  the  plasma  wire  acts  as  If  spring  loaded  by  the 
trapped  compressed  magnetic  field  and  it  expands,  reducing  its  density.  If 
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the  density  would  have  been  high  In  the  first  place  then  most  of  the  described 
action  would  have  had  a  lesser  effect.  We  thus  recotmtend  a  way  of  reducing 
early  plasma  expansion  to  a  minimum.  The  basic  Idea  Is  the  elimination  of 
the  early  field  emission  process.  This  can  be  achieved  by  making  the  wires 
be  at  a  positive  potential  compared  to  the  return  current  rods.  While,  one 
possibility  to  achieve  a  change  of  polarity  is  reversal  of  generator  polarity 
we  do  not  recommend  this  because  of  possible  power  reduction.  The  other  ways 
are  either  using  convoluted  geometry  (Fig.  1)  or  making  a  reversed  arrange¬ 
ment  of  geometry  Inside  the  cathode  (Fig.  2). 


Fig.  1.  Convoluted  cathode  anode. 


Fig.  2.  Wire  array  Inside 
cathode. 


In  these  experiments  the  main  feature  to  be  tested  Is  the  achievement 
of  denser  plasmas.  The  physics  Included  In  the  next  section  should  also  be 
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considered  before  performing  the  suggested  reversed  polarity  experiments  and 
It  Is  desirable  to  check  them  first. 


(2)  Radial  Collapse  of  Wire  Arrays. 

The  present  way  of  thinking  about  the  radial  motion  of  the 
wires  onto  the  diode  axis  assumes  that  the  attraction  between  the  wires 
pulls  them  Inwards.  One  then  arrives  at  a  very  simple  conclusion  that  since 
current  increases  linearly,  I-t,  and  attraction  force  increases  as~I2,  the 
initial  motion  Is  given  by  the  acceleration  equation 


k  t2 


so  that 

r  ~  r  O  -  (|~n  . 
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This  result  will  change  drastically  if  we  look  a  little  bit  closer  at 
the  experimental  arrangement  along  the  wire  array  (.Fig.  3)  taking  as  an 
example  6  wires. 


In  the  typical  experiment  the  distance  between  wires  dg  is  much 
larger  than  d^  (dg/d^  ~  5  to  10).  The  current  outside  does  not  return 
unlformally  through  a  cylinder  generating  no  B  fields  inside.  It  goes 
through  the  rods  and  because  of  skin  effects  runs  at  the  surface  of  these 
rods.  The  single  rod  near  a  wire  generates  a  magnetic  field  larger  than  all 
the  other  wires.  Since  the  currents  In  wire-rods  are  antiparallel,  they 
repel  each  other  and  the  wire  Is  pinched  Inside.  All  the  old  calculations 
have  missed  this  point  and  new  calculations  should  be  performed  to  see  how 
much  of  an  effect  on  assembly  time  will  be  predicted.  In  sunmary:  one 
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Fig.  3.  Wires  and  Return  Current  Rods 


should  say  that  at  early  times  (dg/dj  »  1)  the  wires  are  pushed  by  the  rods 
inside  but  later  In  time  (dg/dj^l)  they  pull  each  other. 

(3)  Stability  of  the  Run  In  Wire  Array. 


We  now  look  at  the  question  of  how  stable  the  wire  array  Is  for 
assembly  on  axis.  This  work  Is  very  different  from  the  one  already  performed 
by  Rostoker  in  a  number  of  aspects.  N.  Rostoker  considered  only  the  attrac¬ 
tive  forces  between  the  wires  and  neglected  the  rods  and  as  a  result  was 
completely  insensitive  to  an  azymuthal  transformation  of  the  whole  wire 
array.  In  this  treatment  one  considers  the  main  effect  to  be  that  of  the 
rods  and  thus  shows  that  any  tiny  transformation  is  Inherently  unstable 
because  of  the  present  shape  of  the  rods.  Actually,  It  Is  difficult  to  con¬ 
ceive  of  a  geometry  that  would  be  more  unstable.  Consider  the  geometry  of 


9 


25 


Fig.  3  and  notice  that  the  current  In  the  rods  runs  very  near  to  Its  tip. 

If  wires  are  not  perfectly  aligned  with  each  other,  than  the  pushing  force 
from  the  rods  will  not  push  them  Into  the  center  and  they  will  only  partially 
collide  there,  overlapping  their  plasmas.  If  they  are  all  aligned  very  well 
with  each  other,  than  the  pushing  force  from  the  rods  will  not  push  them 
Into  the  center  and  they  will  only  partially  collide  there,  overlapping  their 
plasmas.  If  they  are  all  aligned  very  well  with  each  other  then  a  tiny  twist 
of  the  wire  array  relative  to  the  rods  will  result  in  a  rotational  momentum 
causing  the  wires  to  miss  each  other  near  axis. 

This  may  explain  the  large  scatter  In  the  data  of  radiation  for  given 
electrical  pulses  on  the  generator  that  have  basically  the  same  pulse  shape 
and  power.  Viewing  the  wire  array  x-rays  from  the  anode  direction  as  a  func¬ 
tion  of  time  (a  few  frames)  near  assembly  may  show  If  indeed  such  imperfect 
radial  motion  exists  In  an  experiment  where  a  relative  twist  was  intentionally 
superimposed.  It  Is  also  recommended  to  check  time  integrated  data  of  shots 
already  taken  and  see  If  low  radiation  output  was  accompanied  by  non  radial 
Implosion. 

The  question  of  how  to  stabilize  such  an  implosion  Is  treated  next. 

Out  of  the  many  possibilities  it  was  found  that  the  cup  arrangement 
(Fig.  4)  was  the  simplest  and  easiest  experimentally.  The  two  sharp  sides 
act  as  stabilizers  as  they  push  the  wire  radially  Inwards. 

C.  Relativistic  Electron  Beams  in  the  Plasma. 

When  the  flow  of  relativistic  electrons  enters  the  imploding 
plasma  It  participates  In  the  structure  of  the  total  current  and  the  energy 
flow  balance  as  well  as  In  the  pressure  balance.  The  origin  of  these  electrons 
could  be  In  the  larger  radius  feed  and  would  be  magnetically  Insulated  and 
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Fig.  4.  Old  and  New  Return  Current  Rods 


driven  into  the  plasma  by  an  Ef  x  BQ  drift.  When  they  enter  the  plasma 
these  radial  fields  drop  by  a  very  large  factor  due  to  the  high  conductivity 
of  the  plasma  and  its  large  radial  extent  (~1  cm).  Outside  of  the  plasma 
fields  of  the  order  of  1  MV/cm  exist  but  inside  it  they  drop  significantly. 
An  upper  limit  to  the  axial  drift  of  such  electrons  due  to  the  Ep  field  may 
thus  be  estimated  by  assuming  Ef»  3x10s  V/ cm  and  find  that: 

E  x  B  E 

<0  *  c-X—S  =  e  _£  0) 


which  is  about  5xl0“*c  for  currents  of  3  MA  at  3  mm  radii. 

Another  drift  that  takes  place  is  that  of  VB  type.  Very  near  the  axis 
where  the  electron  flow  encloses  to  less  than  the  Alfven  current  (50  kA  for 
1  MeV)  electrons  propagate  by  sinusoidal  orbits  with  fast  average  speeds 
(v~  c)  and  thus  make  It  to  the  anode  losing  almost  no  energy.  The 
fraction  of  energy  observed  using  stacked  films  for  hard  x-ray  observation 
as  reported  by  the  Maxwell  Laboratory  group  (private  communication  with 
Wayne  Clark)  agrees  well  with  one  Alfven  current  of  relativistic  electrons 
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arriving  at  the  anode.  As  one  follows  electron  orbits  at  larger  radii  the 
electrons  drift  via  gradients  In  the  magnetic  field  structure.  The  drift 
is  towards  the  cathode  for  increasing  BQ  fields  and  towards  the  anode  as  the 
field  decreases.  The  value  of  the  drift  velocity  in  (see  Figure  5). 
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where  m  is  the  relativistic  mass  of  the  electron  and  Vj_is  its  perpendicular 
velocity.  In  addition,  curvature  B  drifts  also  are  directed  towards  the 
anode  and  have  a  value  comparable  to  that  of  Eq.  (1)  due  to  azymuthal 
velocities  of  the  relativistic  electrons  thus  one  should  use  twice  the  value 
of  Eq.  (1).  The  magnetic  field  falls  off  slower  than  1/r  on  Its  shoulder  and 
thus  an  overestimate  at  <VZ>  will  be  given  using  an  1/r  dependence.  It  is  then 
easy  to  obtain.  _ _ 


<vz>  < 
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The  electrons  thus  drift  very  slowly  with  tight  orbits  and  the  deposition  of 

their  energy  a-ong  the  axial  direction  Is  enhanced  by  the  ratio  of  =  R 

which  is  larger  than  - .  For  1  MeV  electron  (0  =  0.95,  5x10^) 

81 A1  fven  AITven 

and  for  I  «  3  MA  then  gives  2  x  102  enhancement  of  the  stopping  power  of  the 
electrons  In  the  plasma  (typically  2  MeV/g/cm2).  While  the  diluted  density 
p~10  g/cm3  of  the  imploding  plasma  would  have  subtracted  over  a  distance 
of  3  cm  about  6  keV  from  the  relativistic  electron,  when  the  enhancement 
factor  Is  Included  all  the  1  MeV  energy  Is  lost  In  the  plasma.  We  thus  con¬ 
clude  that  imploding  plasmas  with  large  currents  form  very  good  traps  for 
the  relativistic  electrons  as  long  as  the  Impedance  of  the  generator  is  less 
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than  a  critical  value  so  that  the  energy  loss  in  the  10  g/cm3  plasma  over 
X,  ■  3  cm  length  exceeds  the  voltage  of  the  electrons  (nearly  that  of  the 
diode) 

p  l  R  S  >  V  , 

or  for  the  case  considered 

0.7  I  >  V  B2  Y  .  (4) 

Approximating  y  a  1  +  2V  by  2V  where  V  is  in  MV,  one  obtains 

0.35  I (MA)  >  V2(MV)  . 

The  statement  is  that  If  the  voltage  of  the  feed  is  about  1  MV  then 
one  needs  about  3  MA.  If,  however,  3  MV  are  used  then  about  30  MA  are  needed. 
Otherwise,  the  beam  can  make  it  all  the  way  through.  That  is  probably  the 
reason  why  AMP,  that  performed  as  a  high  voltage,  low  current  generator  showed 
a  lot  of  bremsstrahlung  x-rays  on  the  anode,  causing  only  a  fraction  of  the 
current  to  be  carried  by  the  slow  plasma  electrons.  Low  impedance  generators 
thus  have  inherent  advantages  of  being  able  to  absorb  the  kinetic  energy  of 
the  beam  back  into  the  plasma.  One  word  of  caution  regarding  the  use  of  the 
above  numbers.  The  stopping  power  of  the  relativistic  electrons  is  for  non- 
lonized  medium  that  gives  the  2  MeV/g/cm2.  When  the  plasma  is  very  hot  the 
classical  stopping  power  changes  and  increases  by  factors  of  2  or  3  (depending 
on  Ionization  state,  temperature,  density,  etc.).  This  is  only  a  theoretical 
prediction  since  no  one  has  yet  measured  that  increase.  In  addition,  one  must 
consider  the  magnetic  effect.  The  numbers  In  Eq.  (.4)  may  thus  be  alleviated 
and  it  may  look  like  I  >  V2. 
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Fig.  5.  Magnetically  insulated  REB  trapped  in  the  plasma 
D.  Impact  on  Generator  Design. 

Since  the  voltage  on  the  upper  part  of  the  diode  in  vacuum  is 
given  by  the  expression  V  a  ^  (L*I)  where  L  is  the  diode  Inductance  in 
vacuum  and  since  this  Inductance  does  not  change  much  in  the  beginning  we 
find  that  V  ■  L  ^  .  The  generator  dictates  that  voltage  by  its  open  circuit 
pulse  VQ  c  and  Its  Impedance  ZQ.  Thus  the  voltage  is  V  *  VQ  c  (t)  "  1  zg  ' 

L1nsulator  *  For  given  vo.c.^*  and  can  56  so1ved  exactly' 

Our  new  concern  now  Is  that  the  resultant  current  and  voltage  satisfy  the 

inequality  of  Eq.  (4)  of  Section  C,  which  Is  putting  a  constraint  on  either 

L  or  V  (t).  Specifically  one  should  lower  L  to  a  minimum  (this  will 
o»  c* 

be  discussed  In  the  next  section  on  magnetic  insulation)  but  more  important, 
one  needs  to  give  V.  .  (t)  a  slow  Increase  in  time.  While  the  exact  treat- 
ment  is  rather  lengthy  the  major  results  can  be  obtained  by  simply  assuming 
V  *  L  £  where  x  is  the  rise  time  and  now  use  Eq.  (4)  of  Section  C  to  obtain 
0.3  I  >  (—)*  or  x  >  1.7  L. I**  which  states  that  for  larger  currents  one  likes 
to  have  longer  rise  times,  more  than  100  ns  for  20  nH  and  10  MA.  This  agrees 
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with  the  previously  unexplained  experimental  observation  that  better  results 
were  obtained  with  long  pulses  Including  the  so  called  high  power  prepulse 
of  BJ  5.  In  conclusion,  the  generator  equations  above  need  to  be  solved  for 
given  V.  _  (t)  and  the  resultant  V(t),  I ( t)  checked  to  see  If  inequality  In 

w  •  V  • 

Eq.  (4)  of  Section  C  Is  satisfied,  otherwise  V.  .  (t)  must  be, changed. 

u  •  c  • 

E.  Problems  In  Magnetic  Insulation. 

The  magnetic  insulation  of  the  vacuum  feed  is  important  since  it 
prevents  electron  flow  from  crossing  the  A-K  gap  at  large  radii.  Unlike 
low  voltage  (<  1  MeV)  generators  which  do  not  turn  on  the  cathode,  the  high 
power  generators  operate  with  all  cathode  surfaces  turned  into  a  plasma.  Two 
major  effects  should  be  looked  at:  the  electron  current  loss  to  the  anode 
and  the  ion  current  loss  to  the  cathode.  Specifically,  the  electron  loss  is 
important  since  it  heats  the  anode  surface  outside  the  plasma  load.  This 
heating  will  cause  gas  desorption  off  the  anode  material  igniting  an  anode 
plasma.  The  time  scales  for  such  a  phenomenon  are  derived.  The  ion  current 
is  far  In  excess  of  the  simple  Child-Langmuir  value  due  to  the  presence  of  the 
electron  radial  flow  similar  to  that  In  a  pinched  beam  diode.  For  beam 
electron  flow  of  2 !2  the  ion  current  is  nearly  twice  the  electron  current 
causing  the  total  current  loss  to  be  about  In.  Thus,  if  a  feed  operates  at 
3  MV  and  has  a  total  current  of  10  MA  with  1.5  MA  of  that  current  (only  15%) 
in  the  electron  beam,  this  current  will  cause  (if  allowed  to  ignite  an  anode 
plasma)  an  ion  current  of  3  MA  leaving  only  5.5  MA  in  the  plasma  load  on  axis. 

It  is  therefore  obvious  that  the  electron  current  in  the  beam  should 


be  minimized  and  the  Ignition  of  anode  plasma  should  be  prevented  at  large 
radii  by  choosing  correct  geometries  and  anode  materials  that  can  withstand 
heating  for  longer  times.  First  we  show  the  basic  mathematical  treatment  of 


the  magnetic  insulation  in  a  radial  feed,  than  we  compute  anode  heating  and 
point  out  material  response  for  anode  plasma  ignition  times. 


Electron  orbits  in  radial  feeds 

We  assume  that  a  constant  voltage 
exists  between  anode  and  cathode 
during  the  time  when  the  plasma 
runs  inwards  and  the  effective 
inductance  change  is  such  that 


which  is  satisfied  by  most  low 
Impedance  high  power  generators. 


Fig.  6.  Typical  electron  orbit  In  vacuum  feed 

The  geometry  considered  and  a  typical  electron  orbit  are  shown  In  Fig. 

6  .  The  electron  drifts  radially  Inwards  by  the  ExB  drift  and  It  drifts 
axially  across  by  the  VB  drift.  These  two  drifts  dictate  a  given  path  for 
the  guiding  center  motion  provided  the  axial  (the  A-K  gap)  and  radial 
(radius  R)  lengths  are  much  larger  than  the  larmour  radius  rL  ■  1  ~B(kG) U —  • 
Since  the  typical  A-K  gaps  are  one  centimeter  and  since  currents  larger  than 
3  MA  will  be  considered,  the  typical  magnetic  fields  at  radii  on  the  order  of 
10  cm  are  larger  than  60  kG.  The  gyroradll  are  of  the  order  of  millimeters, 
much  smaller  than  A-K  gaps.  The  equations  for  the  axial  and  radial  velocities 
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are: 


2 
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Here  we  have  assumed  that  the  E  field  is  a  constant  in  the  gap,  E  *  g— » and 
B~r-1.  Also  y  is  the  relativistic  mass  ratio  for  the  electron  which  is  a 
function  of  its  position  and  VQ  and  I  are  the  voltage  and  current  in  MV  and 
MA.  Note  that  the  radial  velocity  is  inversely  proportional  to  the  impedance. 


The  equation  of  the  orbits  is  derived  through 


c 


(3) 


and  the  result  is 

dr  -  r 
dZ  “  1 

Y“  “ 


5a 
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Since  v~l  +  2  V^-l  +  2VQ  |  ,  we  may  separate  the  equation 


dZ 


4V 
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and  Integrate 


(4) 
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Here,  the  particle  started  at  a  distance  from  the  cathode  at  a  radius  r^ 
and  will  reach  a  distance  Z  at  a  radius  r  given  by 
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(5) 


r,  7  V  V 

-i,Z  <)+  <>z)  (,♦  °  z) 


The  question  of  finding  the  radius  where  the  particle  crosses  the  gap  D  can 
now  be  answered: 


r  = 


7  -2- 

h  0  +  d  zo 

D  1  +  VQ 


(6) 


The  result  shows  that  large  gaps  D  are  very  beneficial  and  also  shows 
how  sensitive  the  result  is  on  the  value  of  Z-|  and  r^  namely  the  specific 
flow  pattern  that  lift  electrons  off  the  cathode  shank  dictates  how  well  they 
can  be  magnetically  insulated  flowing  radially  inwards.  There  is  a  definite 
advantage  of  positive  polarity  operation  since  it  eliminates  shank  structure 
and  thus  minimizes  Z-j.  For  simplicity  we  assume  that  electrons  lift  off 
by  the  electric  field  so  that  near  the  cathode  all  the  velocity  gained  turns 
parallel  to  it  at  a  distance  2  Z-|  (see  Fig.  7). 


Figure  7 
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Near  y~l,  one  may  write  the  following  expressions 

p2  +  p2 

(for  relativity)  y2  *  1  +  — - ~  »1  +  B2  +  g2 

^oc)2  "  2 


or  y2  ~  1  +  82  at  2 

V 

(for  energy)  y  =  1  +  2  VMy(Z1)  ^  1  +  4  Z]  ^ 

eB  22 

(for  the  radial  momentum  due  to  B)  8  « - 

r  me2 

After  some  algebra  and  using  B  ~  — ---t— )  ,  these  reduce  to 

rl 


(7) 

(8) 
(9) 


ll  -  Vl  1-5  x  10~4 
»'1  D  1 2  ( MA) 


(10) 


Inserting  (10)  into  (6)  one  finds 


r 


Vo  1.5  x  10~4 
V0  I2  (MA) 


1.5  x  10~4 
P  ( MA) 


(ID 


The  dependence  on  VQ  Is  weak,  the  very  large  currents  help  but  large  radii 
are  the  biggest  problem.  One  must  design  around  it  to  try  to  minimize  it. 

For  the  typical  values  I  -  3  MA,  r}  »  100  cm,  and  D  »  2  cm,  one  obtains 

r  »  4  cm  (which  can  be  made  1  cm  If  D  =  4  cm).  For  I  *  10  MA  and  all  of  the 

above  the  same,  r  =  0.4  cm  which  brings  us  back  to  trapping  the  electrons  in 

the  Imploding  plasma.  Changing  radii  to  r1  *  200  cm  causes  a  problem  again 
since  r  *  4  cm. 
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F.  Igniting  the  Anode  Plasma. 

Assuming  now  that  a  finite  fraction,  Ifc,  of  the  total  current  flows 

as  a  beam  Igniting  an  area  up  to  radii  of  ~ 5  cm,  we  find  that  an  area  A  of 

about  75  cm2  is  heated.  If  one  uses  a  given  anode  material  the  heating  will 
Ignite  a  plasma  (independent  of  electron  voltage)  on  a  time  scale,  t,  given 
by  using  the  concept  of  areal  velocity  developed  for  the  study  of  pinched 
beam  diodes, 

\  r  1  =  1  •  cia) 

In  order  to  increase  t  (so  no  ion  current  flows  during  it)  one  has  to 

use  materials  with  small  VA  such  as  carbon  or  coatings  of  teflon  or  polyethelyne 

that  all  have  high  specific  heat,  causing  longer  times  for  needed  heating. 

If,  however,  one  uses  stainless  steel,  its  low  specific  heat  increases  by 
2 

VA  *  10  lOO^OTns  which  for  lb  =  300  kA,  A  *  75  cm2  gives  t  =  2.5  ns  while 

for  carbon  it  will  be  about  50  ns.  If  ion  beams  are  allowed  they  are  I^s*2I@ 

for  ~lfl  structures  and  thus  about  1  MA  will  be  lost  In  the  above  example. 

The  above  examples  obviously  point  out  the  Importance  of  structure  and 
anode  materials  and  their  effects  on  the  existence  of  mechanisms  for  current 
flow  in  the  diode. 
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